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Abstract 

The breaking of time- reversal symmetry by ferromagnetism is predicted to yield profound changes 
to the electronic surface states of a topological insulator. Here, we report on a concerted set of 
structural, magnetic, electrical and spectroscopic measurements of Mn-Bi2Se3 thin films wherein 
photoemission and x-ray magnetic circular dichroism studies have recently shown surface ferro- 
magnetism in the temperature range 15 K < T < 100 K, accompanied by a suppressed density 
of surface states at the Dirac point. Secondary ion mass spectroscopy and scanning tunneling 
microscopy reveal an inhomogeneous distribution of Mn atoms, with a tendency to segregate to- 
wards the sample surface. Magnetometry and anisotropic magnetoresistance measurements are 
insensitive to the high temperature ferromagnetism seen in surface studies, revealing instead a low 
temperature ferromagnetic phase at T < 5 K. The absence of both a magneto-optical Kerr effect 
and anomalous Hall effect suggests that this low temperature ferromagnetism is unlikely to be a 
homogeneous bulk phase but likely originates in nanoscale near-surface regions of the bulk where 
magnetic atoms segregate during sample growth. Although the samples are not ideal, with both 
bulk and surface contributions to electron transport, we measure a magnetoconductance whose be- 
havior is qualitatively consistent with predictions that the opening of a gap in the Dirac spectrum 
drives quantum corrections to the conductance in topological insulators from the symplectic to the 
orthogonal class. 
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I. INTRODUCTION 



The metallic Dirac cone surface states of a topological insulator (TI) are expected to be 
protected against small perturbations by time-reversal symmetry^. However, these surface 
states can be dramatically modified by breaking time-reversal symmetry through magnetic 
doping^. This is predicted to open a finite energy gap at the Dirac pointP, an important 
consideration for potential device applications of TIs. Further, the interplay between mag- 
netism and TI surface states is predicted to yield novel phenomena of fundamental interest 
such as a topological magneto-electric effectP, a quantized anomalous Hall effectP, and the in- 
duction of a magnetic monopole 7 . These predictions have already triggered an active search 
for materials that provide experimental realizations of such phenomena. Angle-resolved 
photoemission spectroscopy (ARPES) studies have shown that magnetic impurities modify 
the surface states^ and may possibly open a surface state gapP. Ferromagnetism has been 
demonstrated in several magnetically doped TI crystals, including bulk crystals of Mn-doped 
BiaTe^and B^Tea^Se™, as well as thin films of V-doped Sb 2 Te^, Cr-doped Sb 2 TeJ™ 
and Cr-doped Bi 2 Se^H The evidence for ferromagnetism in these previous studies is pri- 
marily based upon magnetometry and magnetotransport, neither of which unfortunately 
provide a clear means of distinguishing between surface and bulk magnetism. However, the 
behavior of the anomalous Hall effect in electrically gated Mn-doped Bi2Te3_ 3/ Se y samples 
has been interpreted as evidence for surface ferromagnetismfm. 

Recently^, a spin resolved ARPES and x-ray magnetic circular magnetic dichroism 
(XMCD) study of Mn-doped Bi 2 Se3 thin films clearly revealed the presence of out-of-plane 
surface ferromagnetism that persists up to relatively high temperature (T ~ 100 K) and 
demonstrated a hedgehog spin texture accompanying a Dirac-metal to gapped-insulator 
transition. The broadening of the ARPES data also suggested an inhomogeneous surface 
magnetization. Here, we report additional studies of these thin films down to low tem- 
peratures (T > 300 inK), focusing on concerted structural, magnetic, electrical, and spec- 
troscopic measurements. Magnetometry and magneto-transport measurements indicate the 
existence of a low temperature ferromagnetic phase with a Curie temperature Tq < 5.4K 
and an in-plane magnetization. The use of multiple complementary structural probes in our 
Mn-Bi 2 Se3 thin films suggests that this observed low temperature ferromagnetic phase is 
not a homogeneous bulk effect but likely occurs in a limited region (~ 5nm) near the sam- 
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pie surface. We note that these macroscopic magnetization and transport measurements 
are likely to miss the purely surface out-of-plane and high Tq ferromagnetism observed in 
spin resolved ARPES measurements and XMCD. We also note that a recent calculation^ 
indicates that strictly surface ferromagnetism can be more robust than ferromagnetism me- 
diated by bulk states, a scenario that is consistent with our observations. As the samples are 
cooled through Tc of the low temperature ferromagnetic phase, the magnetoconductance of 
the samples goes through a distinctive transition from symplectic (or weak anti localization) 
to orthogonal (or weak localization) quantum corrections to transport. Since ARPES mea- 
surements show that the Fermi level in these samples is located well in the bulk conduction 
band states, the magnetoconductance data probe a complicated interplay between states 
that include the pure TI surface, confined two dimensional (2D) states created by band 
bending and three dimensional (3D) bulk states. Despite this complex situation, the onset 
of time reversal symmetry breaking ferromagnetism qualitatively yields the magnetoconduc- 
tance behavior predicted for 2D Dirac cone TI surface transport^. We note that a recent 
study 19 of ultrathin (3 QL) samples of Cr-doped Bi 2 Se3 showed a similar cross-over from 
weak anti-localization to weak localization, but in a regime where the samples are paramag- 
netic and where a surface state gap is already developed at a temperature (~ 100 K) much 
higher than that implied by transport measurements (< 10 K). Given the development of a 
gap due to the hybridization of the upper and lower surfaces in these ultra thin samples, it is 
difficult to unambiguously disentangle this effect from the opening of an additional gap from 
possible breaking of time reversal symmetry. Nonetheless, it is quite possible that surface 
ferromagnetism also exists in those samples and could be detected using XMCD. 

II. THIN FILM SYNTHESIS AND CHARACTERIZATION 

The Mn-Bi2Se3 films were synthesized by MBE using high purity elemental (5N) Mn, Bi, 
and Se sources. After thermal desorption of the native oxide on the epi-ready GaAs 111 A 
substrate, a thin ZnSe buffer layer (~ 8 nm) was first deposited at ~ 300°C using a Se/Zn 
beam equivalent pressure (BEP) ratio of ~ 2.3 to yield a smooth surface. A Mn-Bi 2 Se3 layer 
(~ 60 nm) was then synthesized immediately on the ZnSe buffer with a Se/Bi BEP ratio of 
~ 15. The Mn doping concentration was controlled by changing the Bi/Mn BEP ratio in a 
range of ~ 8 — 24. Reflection high energy electron diffraction (RHEED) was used to monitor 
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the entire growth process. Figure 1 (a) shows a typical atomic force microscopy (AFM) 
image of a Mn-Bi 2 Se3 thin film (BEP ratio Bi/Mn = 12.5) with a root mean square surface 
roughness ~ 0.5 nm. The triangular islands at the surface indicate screw dislocations. X- 
ray diffraction (XRD) shows that Mn-Bi2Se3 has the same rhombohedral layered crystal 
structure as un doped Bi2Se3, but with strained lattice parameters. The observation of only 
(003) family of reflections indicates that the growth direction is oriented along the c axis 
with no secondary phases such as MnBi or MnSe (Fig. 1 (b)). In highly doped samples, the 
smaller interplanar distance between (006) planes indicates that the Mn atoms substitute 
the Bi atoms due to the smaller atomic radius of Mn (0.140 nm) relative to Bi (0.160 nm); 
while the larger interplanar distance between (003) planes indicates that the Mn atoms also 
intercalate between quintuple layers. As the Mn doping concentration increases, the quality 
of the thin films degrades and the crystal structure shifts towards the BiSe phase. 

Inhomogeneous phase separation into nanoscale regions of ferromagnetism is well known 
in magnetically doped semiconductors^. To examine this possibility, we measured the sec- 
ondary ion mass spectrometry (SIMS) depth profile in several samples to reveal the Mn 
concentration distribution along the growth direction. Figure 1 (c) shows a depth profile 
from the Bi/Mn = 12.5 sample. The Mn doping concentration (~ 20% relative error) was 
determined using a sample with bulk Mn concentration of x = 0.043 (Mn r (Bi2Se3)i_ ;c ) as 
a standard, the composition of which was obtained from Rutherford backscattering spec- 
troscopy. These measurements show that the surface Mn concentration is almost an order of 
magnitude higher than that of the bulk. X-ray photoelectron spectroscopy (XPS) measure- 
ments provide further confirmation of excess Mn at the surface: the averaged Mn concentra- 
tion from the top ~ 5 nm of the Bi/Mn = 12.5 sample is x = 0.085, about eight times higher 
than that of the bulk (x = 0.010). To study the morphology of the excess Mn at the surface, 
we carried out scanning tunneling microscopy (STM) measurements on samples that were 
capped in situ with Se immediately after growth. Careful de-capping of these samples in an 
ultrahigh vacuum chamber at 200 °C, lower than the growth temperature of 300 °C, allows 
the acquisition of STM data that reveal substitutional Mn below the surface, individual 
Mn adatoms on the surface, as well as numerous small regions with a high density of clus- 
tered Mn atoms (Fig. 1 (d)). The clusters have typical dimensions of tens of nanometers, 
are generally located near atomic step edges (indicated by dark lines), and are distributed 
across the entire surface. Cross-sectional high-resolution transmission electron microscopy 
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(HRTEM) provides further detailed insights into the bulk structural characteristics of the 
Mn-Bi 2 Se 3 thin films. Figures 1 (e) and 1 (f) show a representative HRTEM image from a 
cross-sectional specimen of the Bi/Mn = 12.5 sample and the corresponding selected area 
electron diffraction pattern, respectively. The clear phase contrast image and diffraction 
pattern indicate that the Mn-Bi2Se3 films are generally high quality single crystal with little 
disorder. Within our TEM resolution limit, there is no clear evidence for the formation of 
Mn clusters in the bulk. The observation of Mn clusters at the surface in STM studies and 
their absence in the bulk in HRTEM studies further confirm that there is excess surface Mn 
in these samples. 

III. MAGNETOMETRY MEASUREMENTS: OBSERVATION OF FERROMAG- 
NETISM 

After understanding the Mn distribution in our Mn-Bi 2 Se 3 thin films, we address the 
possibility of time reversal symmetry breaking ferromagnetic order amongst the Mn impu- 
rities. As we mentioned in the introduction, XMCD measurements have indeed revealed 
the presence of surface ferromagnetic order with out-of-plane magnetization up to at least 
T = 100 K. We have now carried out magnetization measurements using a superconducting 
quantum interference device (SQUID) magnetometer in a Quantum Design Magnetic Prop- 
erty Measurement System with a base temperature at 1.8 K. Figure 2 (a) shows temperature 
dependent magnetization of the Bi/Mn = 12.5 samples in field along both the (1100) direc- 
tion and the c axis. The temperature dependent magnetization curves were obtained using 
a 50 Oe measuring field after field cool with H = 9 kOe. Fitting to the high-temperature 
data allows for the determination of a temperature-independent diamagnetic signal of the 
substrate, the constant offset of which is removed from the sample data. Fitting a line to 
the steepest slope of the curve allows us to extrapolate a ferromagnetic Curie temperature 
T c ~ 5.4 K (indicated by the arrow). Based on the Mn distribution in this thin film, we es- 
timate ~ 3.45/iB per Mn atom by assuming that the ferromagnetism originates from the top 
5 nm near the surface with 8.5 a.t.% Mn based on our XPS measurement. Figures 2 (b) and 
2 (c) plot magnetization vs. field curves of the Bi/Mn = 12.5 sample at different tempera- 
tures in fields along the (1100) direction and the c axis, respectively. The diamagnetic signal 
due to the substrate is subtracted off by fitting to the high-field data and correspondingly 
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removing the diamagnetic background. The hysteresis above T c (Fig. 2 (c)) is attributed 
to defects and/or impurities in the GaAs substrate based on careful control measurements 
of several bare GaAs substrates. The easy axis is found to be in plane and close to the 
(1100) direction. Ferromagnetic ordering with a Tc ~ 3.1 K is also observed in the Bi/Mn 
= 23.6 sample (Fig. 2d - 2f). In the two lightly doped samples (Bi/Mn = 23.6 and 12.5), 
Tc increases with Mn doping, a trend that we have noticed in other samples grown in this 
doping range (data not shown). However, we also find that for high Mn doping, the quality 
of the crystal degrades and the crystal structure shifts towards BiSe instead of Bi 2 Se 3 . A 
higher Mn concentration also increases the carrier concentration significantly (up to ~ 10 20 
cm -3 ). Both effects seem to result in a lower Tc- 

IV. TRANSPORT MEASUREMENTS: OBSERVATION OF FERROMAG- 
NETISM AND CROSS-OVER FROM WEAK ANTI-LOCALIZATION TO WEAK 
LOCALIZATION 

The ferromagnetic ordering in Mn-Bi 2 Se3 thin films is further confirmed by magneto- 
transport measurements. To study electrical transport, Hall bar devices (1300 /im x 400/zm) 
were fabricated from thin films using photolithography and wet etching. Transport measure- 
ments were carried out using a standard lock-in technique in an Oxford Heliox He-3 cryostat 
with a base temperature at 0.5 K and field up to 6 T as well as an Oxford Triton He-3 vector 
magnet system with a base temperature at 0.3 K wherein a 1 T field could be arbitrarily 
oriented over the entire sphere. For clarity, we focus on the sample with Bi/Mn = 12.5, 
noting that similar behavior is also observed in other samples. The measurement geometry 
is shown in Fig. 3 (a): the Hall bar device lies in the xy plane with the excitation current 
along the x axis ((1100) direction). The orientation of the magnetic field H is defined by the 
polar angle 9 and the azimuthal angle (p. Figure 3 (b) shows the magnetoconductivity (MC) 
Aa xx = a xx (H) — a xx (0) in units of e 2 /h at T = 0.3 K when the magnetic field is swept at 
different fixed angles in the xz plane. The arrows indicate the field sweep directions. The 
longitudinal conductivity a xx is obtained from the longitudinal and transverse resistivity 
measurements using the 2D tensor: a xx = ( p 2 p + p 2 y Note that hysteresis is observed in the 
MC for all field orientations, but is most pronounced when the magnetic field is at a slight 
angle from the z axis. This angular dependence of the hysteresis in the MC is a strong indi- 
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cation of ferromagnetic ordering amongst the magnetic impurities, with a clear anisotropy 
that favors a magnetization oriented close to the x axis. To further identify the critical tem- 
perature T c of ferromagnetism, we carried out temperature dependent measurements of the 
MC with the field oriented at the angle that showed the most pronounced hysteresis (6 = 5° 
and if = 0). Figure 3 (c) shows that the magnitude of the hysteresis and the magnetization 
switching field both decrease with increasing temperature, vanishing at T c ~ 5.5 K, which is 
consistent with the T c obtained from SQUID measurement (Fig. 2 (a)). Similar hysteretic 
behaviors were also observed in other samples (Fig. (4)). 

The presence of ferromagnetism is also manifest through another classic signature in 
transport: anisotropic magnetoresistance (AMR). In these measurements, we carried out 
angular sweeps of a magnetic field with fixed magnitude in the xy plane. In Fig. 3 (d), 
an in-plane magnetic field was applied with a constant magnitude of 0.5 T to saturate 
the magnetization while its orientation was continuously swept through 360°. The angu- 
lar dependence of the magneto-resistance can be interpreted using the standard heuristic 
expression for AMRP, 

E x = jp± + j(p\\ - p±) cos 2 if, (1) 

where E x is the longitudinal electric field within a single domain ferromagnetic film, j is the 
current density, <p is the angle between the magnetization and current density, and p\\ and 
p_i_ are the resistivities for current parallel and perpendicular to the magnetization. The fits 
(lines) using equation (1) describe our data (squares) very well. The minimum resistance 
occurs at ip = and 180°, indicating that p\\ < p±. As the temperature increases, the 
angular dependence of the in-plane MR becomes weaker and finally isotropic above ~ 10 
K. The temperature dependence of AMR is qualitatively consistent with that of the MR 
hysteresis. The slightly higher temperature at which the AMR disappears is attributed to 
contributions from MR anisotropy which is present even in undoped B^SegpZl. 

The observation of ferromagnetism using magnetometry and magnetoresistance would 
usually imply the presence of an anomalous Hall effect (AHE). Surprisingly, we do not 
observe any obvious signatures of an AHE in any of the samples (Fig. 5). Furthermore, 
magneto-optical Kerr effect (MOKE) measurements using light of wavelength 690 nm do 
not reveal any observable signal (down to a precision of 50 pr&d) at temperatures as low 
as 2 K. From the absence of an AHE and MOKE signal, it is at first tempting to attribute 
the ferromagnetism in these samples to nanoscale phase-separated Mn-rich regions located 
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on and near the surface, which is supported by the surface and structural characterization. 
While the presence of these clusters would lead to spin-dependent scattering contributions to 
the MC, they would be unlikely to contribute to an AHE or to a measurable MOKE signal. 
We are aware of one other hybrid micro-composite ferromagnet - Cdi_ a; Mn 3 ;GeAs2:MnAs 
- where AHE is absent despite the clear observation of ferromagnetism^. However, on- 
going measurements of ferromagnetic resonance (FMR) in our Mn-Bi2Se3 samples seem to 
suggest otherwise: the angular dependence of FMR shows the standard behavior seen in 
homogeneous ferromagnets^l. We do not currently have a definitive explanation for all these 
incongruous results, but a reasonable scenario is the following. The out-of-plane surface 
ferromagnet ism seen by XMCD at high temperatures is clearly undetectable using SQUID 
magnetometry, MOKE and AHE. As the sample temperature is lowered below ~ 5 K, the 
Mn-rich bulk region within ~ 5 nm of the surface undergoes another ferromagnetic phase 
transition that is seen by SQUID and by FMR, as predicted by mean field theory^. This 
ferromagnetic order is also too weak to be detected by MOKE, perhaps because the photon 
wavelength used (1.8 ev) is far above the band gap of the bulk semiconductor (0.3 ev); 
magneto-optical effects in semiconductors are well known to be resonantly enhanced near 
optical transitions and can be much weaker under non-resonant conditions. The absence 
of the AHE could simply arise because it is masked by the ordinary Hall effect from the 
parallel bulk channel which dominates the transport. 

Having provided evidence for ferromagnetism in these Mn-Bi 2 Se3 thin films, we now 
discuss the influence of this ferromagnetism on the MC in a magnetic field applied perpen- 
dicular to the sample plane (along the c axis). Figures 6a - d show the longitudinal MC 
at various temperatures for four Mn-Bi2Se3 thin film devices. If we focus on their behavior 
at low fields, a generic pattern emerges in all the samples: at high temperatures, all the 
samples show a negative MC which crosses over to a positive MC as the temperature is 
lowered. The cross-over coincides with the onset of ferromagnetic hysteresis in the MC. In 
contrast, the MC behavior of un doped Bi 2 Se3 is always dominated by a negative MC that 
has a characteristic cusp-like form at low fields, originating from weak anti-localization, and 
a (classical) parabolic or linear form at high fieldsP322 This contrast is apparent in Fig. 6e, 
where the MC response is depicted for samples of various Mn concentrations. 

We now discuss an interpretation of these data based upon diagrammatic calculations^^ 
for quantum corrections to transport in Bi 2 Se 3 . In samples such as the ones studied here, the 
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bulk states dominate the conductivity because of the large Fermi energy: a simple estimate 
using the surface state energy dispersion and a Fermi energy 300 meV above the Dirac point 
shows that the surface carrier density (~ 1.3 x 10 13 cm -2 ) is about an order of magnitude 
smaller than that of the bulk in our samples. Thus, an important question to address is 
the coexisting surface and bulk conduction in our samples. The classical contribution to the 
MC from bulk channels is well known to result in a parabolic positive magnetoresistance or 
a negative MC. According to diagrammatic calculations^, the quantum corrections to the 
conductivity of the lowest 2D bulk quantum well states in Bi 2 Se 3 are purely in the orthogonal 
class and can only result in weak localization (positive MC). In contrast, similar calculations 
for the Dirac surface state show that the quantum corrections to the conductivity are in 
the symplectic class and result in weak anti-localization (negative MC). The calculations 
also show that the opening of a gap at the Dirac point leads to a change in the quantum 
corrections away from the symplectic class, causing a transition from negative MC to positive 
MC with the onset of ferromagnetisnP^. This transition is tunable by the relative size of 
the gap A and the Fermi energy Ep measured from the Dirac point. In the limit of small 
A/2Ep, the MC is dominated by weak anti-localization and is negative. In the limit of large 
A/2Ep, the MC mainly originates in weak localization and is positive. At intermediate gap 
sizes, at low fields a relatively flat classical MC could be expected. While the large number 
of free parameters in the theoretical expressions of ref. 18 precludes meaningful fitting of our 
experimental data, it is still useful to compare the measured MC with simulations obtained 
from the theory. The MC is given byP^ 



where «o and «i are pre-factors for weak localization and weak anti-localization respectively, 
\^(x) is the digamma function, Is is the magnetic length, 1$ is the phase coherence length, 
and /o, h are corrections to l^. Using reasonable parameters (/^ = 300 nm and magnetic 
scattering length l m = 1000 nm) for the simulation, the MC in Fig. 6f shows a cross-over from 
weak anti-localization to weak localization as we adjust A/2Ep from (gapless surface state) 
to 0.99 (fully gapped surface state). This interpretation ascribes all weak anti-localization 
signatures in ungapped topological insulator thin films to surface state transport, regardless 
of the presence of a parallel bulk transport channel. We caution that our films are not thin 
enough to exclude the presence of quantum corrections to the MC from 3D bulk transport 
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as well. These quantum corrections could very well result in a negative MC depending on 
the specific physical parameters^, but we are unaware of a detailed calculation for Bi 2 Se 3 . 

Applying this picture to our experimental observations leads to a consistent explanation. 
The observation of a negative MC at high temperatures (T > 6K) or at high fields is most 
likely due to classical contributions from 2D/3D bulk states. The observed negative MC in 
the Mn-Bi2Se3 samples is not characterized by the typical "sharp cusp" seen in weak anti- 
localization of Bi 2 Se3 samples, though some sign of it can be seen in the magnetoresistance 
above the bulk Tc (Fig. 3 (c)). Since ARPES and XMCD clearly show surface ferromagnetic 
order and a gap at high temperatures, the high temperature curves are likely to correspond 
to the low to intermediate A/2Ep regime, similar to the relatively flat 0.2 ratio curve in 
Figure 6f. As the temperature is lowered (T < 6 K), besides the presence of the out- 
of-plane high temperature surface ferromagnetic phase, a limited near-surface bulk region 
also undergoes a transition to the ferromagnetic phase which is largely in-plane but also 
has a small out-of-plane component (Fig. 2). This in-plane ferromagnetism would not be 
expected to have much effect on the gap except that the coercive field is small enough 
that the application of the out-of-plane external measurement field is enough to rotate its 
orientation out-of-plane. At a field (^ 0.2T) near where the hysteresis loop closes, indicating 
the moments are uniformly pointing out-of-plane, there is a transition to a region of strong 
positive MC consistent with a larger opening of the gap and an increased A/2Ep ratio 
pushing the system into weak localization.^ This interpretation is also consistent with the 
observed absence of MC at low perpendicular fields: in this regime, the magnetization is still 
mainly in-plane and does not contribute much to the further opening of the surface state 
gap, while the disordered out-of-plane moments that do exist could be expected to suppress 
quantum corrections by introducing random phases in the backscattered loops. However, 
once the perpendicular field is large enough, the magnetization begins to be dominated by 
a uniform out-of-plane component, resulting in the positive MC before eventually turning 
over again due to the classical negative MC at high fields (see for example Figs. 6 (a) and 
(b)). The dependence on the Mn doping can also be understood. For the lightly doped 
samples (Bi/Mn ^ 12.5), the carrier concentration is about 4 x 10 19 cm -3 . The increase 
in A with a relatively low Ep gives a larger A/2Ep, hence the positive MC gets enhanced. 
As we further increase the Mn concentration (Bi/Mn ^ 10.3), any further increase in A is 
offset by a considerable rise in the Fermi energy as additional Mn n-dopes the material (the 
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carrier concentration increases to n ~ 2 x 10 20 cm 3 ). As a result, the change in positive 
MC is relatively small compared to the lightly- doped samples. 

V. ARPES MEASUREMENTS 

As we mentioned in the introduction, a detailed study and analysis of ARPES data 
(including spin-resolved measurements) has been reported elsewhere^ and discusses the 
complex spin configuration of the surface states and its time reversal breaking nature as a 
result of magnetic doping. For the sake of completeness and to discuss the relationship be- 
tween the magnetometry, magneto-transport and spin- resolved ARPES, we briefly describe 
the salient features of these surface-sensitive ARPES measurements. These experiments 
were carried out using Se-passivated Mn-Bi 2 Se3 samples (Bi/Mn = 15) which were prepared 
identically to those used in STM experiments. Again, the capped thin films were heated 
up carefully under ultrahigh vacuum (< 1 x 10~ 9 Torr) to remove the Se capping layer. 
ARPES measurements were performed at 20 K on beamline 10.0.1 at the Advanced Light 
Source in Lawrence Berkeley National Laboratory. Typical energy and momentum resolu- 
tion is better than 10 meV and 1% of surface Brillouin zone, respectively. ARPES core 
level spectroscopy data (Fig. 7a) show that only the Se core level was observed before the 
decapping process (blue curve), whereas both Se and Bi peaks appeared after the decapping 
(red curve). Figures 7 (b) and 7 (c) allow a clear comparison between ARPES measured 
dispersion mappings of undoped Bi2Se3 and Mn-Bi2Se3 thin films: while a bright and intact 
Dirac point is observed for the undoped Bi 2 Se3 sample, the energy spectra weight of the 
Dirac node of the Mn-Bi 2 Se3 samples is strongly suppressed. This comparison indicates 
that when time reversal symmetry is broken by introducing magnetic impurities into the 
system, the degeneracy of the electronic states at the Dirac node is lifted. Nevertheless, the 
suppressed intensity at the Dirac point does not go to zero. The absence of a fully developed 
gap may be attributed to the inhomogeneous Mn concentration on the surface, causing a 
spatially dependent magnetic field from the Mn atoms or surface clusters. 

Note that these surface-sensitive ARPES measurements were carried out at temperature 
T = 20 K, which is higher than the Curie temperature observed by magnetization and trans- 
port measurements. However, since the topological surface states are reported to localize 
mainly within the top quintuple layer (also within 1 nm)P™ the electronic structure of the 
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surface states is most relevant to the magnetic properties of the surface. The electronic 
structure within 1 nm of the terminated surface has multiple components that comes with 
spin physics: spin polarized topological surface states whose texture is modified by Mn- 
moments and bulk bands that are bent near the surface which also develop spin polarization 
through Rashba-type effect near the surface. In addition, surface magnetic properties can 
significantly differ from the bulk due to the strongly inhomogeneous Mn concentration pro- 
file along the out-of-plane z direction and the possible surface magnetic ordering metallically 
mediated by topological surface conduction electrons with high mobility. The spin-resolved 
ARPES and XMCD measurements show an out-of-plane magnetic component on the very 
surface (possibly within 1 nm) region of the films under UHV conditions. Therefore, it is 
probably necessary to consider the effects and influence of both the in-plane and the out- 
of-plane magnetic moment components as well as bulk moments of these Mn-Bi2Se3 films 
to understand the full details of magnetism probed in transport and magnetometry. The 
Dirac point spectral weight suppression seen in Fig. 7 (c) originates from multiple physical 
mechanisms: first, the out-of-plane magnetic component opens up a time reversal symmetry 
breaking gap at the Dirac point, which leads to Dirac point spectral weight suppression un- 
der the inhomogeneous Mn distribution on the film surface; second, the in-plane magnetic 
component also leads to spectral weight suppression (without opening a gap at the Dirac 
point). An in-plane magnetization along a fixed momentum-space direction shifts the Dirac 
point away from the T point along the magnetization direction. At T = 20 K, the in-plane 
magnetic moment component on the surface is randomly distributed rather than aligned 
along a fixed direction. Such random in-plane magnetic component consequently leads to a 
substantial broadening of the surface states along the momentum axis. It has been shown 
by several ARPES experiment j 32 l 33 l that substantial momentum broadening of a Dirac-like 
band leads to a spectral weight suppression (a "gap" -like feature) at the Dirac point. Ad- 
ditionally, such observations are also understood by STM measurements^ in the real-space 
point of view by the momentum fluctuation of the surface Dirac fermion in real space. Our 
low temperature transport measurements on Mn-Bi2Se3 thin films at T < 6 K are likely 
related to the combined effect and influence of the in-plane and out-of-plane (purely surface 
origin) magnetic component (the momentum broadening of the surface states) on the shown 
ARPES measurements. A complete picture will require a detailed understanding of the 
interplay between the bulk moments, surface moments and the role of additional Rashba 
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bands near the surface. Due to the inherent surface-sensitive nature of ARPES, currently it 
is not possible to probe bulk magnetism with ARPES. 
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FIG. 1: (Color online) (a) AFM image of a Mn-Bi 2 Se 3 thin film, (b) XRD from four Mn-Bi 2 Se 3 thin 
films, (c) SIMS depth profile showing a higher Mn concentration near the surface, (d) STM 
topographic image (150 x 150 nm 2 , 400mV, 40pA) of the surface of a Mn-Bi2Se3 thin film after 
desorbing the Se capping layer. Substitutional Mn (triangle) clusters around the atomic step edges 
(dark line). Mn diffused onto the surface as adatoms is indicated by the circle and the screw 
dislocation is indicated by the arrow. The insets zoom in on a substitutional Mn atom and a 
Mn adatom (scale bar: 1 nm). (e) Cross-sectional HRTEM image of a Mn-Bi2Se3 thin film, (f) 
Selected area electron diffraction pattern from the same Mn-Bi2Se3 thin film in panel e. 
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FIG. 2: (Color online) (a) Temperature dependent magnetization of the Bi/Mn = 12.5 sample 
in field along (1100) direction and the c axis. The arrow indicates the Curie temperature, (b) 
Magnetization vs. field plots of the Bi/Mn = 12.5 sample in field along (1100) direction at different 
temperatures, (c) Magnetization vs. field plots of the Bi/Mn = 12.5 sample in field along the c 
axis at different temperatures, (d) Temperature dependent magnetization of the Bi/Mn = 23.6 
sample in field along (1100) direction and the c axis. The arrow indicates the Curie temperature, 
(e) Magnetization vs. field plots of the Bi/Mn = 23.6 sample in field along (1100) direction at 
different temperatures, (f) Magnetization vs. field plots of the Bi/Mn = 23.6 sample in field along 
the c axis at different temperatures. 
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FIG. 3: (Color online) (a) Measurement geometry: the Mn-Bi2Se3 Hall bar lies in xy plane with 
the excitation current flowing along x axis. The orientation of the field H is denoted by the polar 
angle 9 and azimuth angle <p. (b) Aa xx versus H plots in different field orientations in xz plane. 
The arrows indicate the field sweep directions and the curves at different field orientations are offset 
for clarity with respect to the 9 = 90° curve, (c) R xx versus H plots in field sweeps (9 = 5° and 
cp = 0) at different temperatures. The hysteresis disappears when weak anti-localization behavior 
starts to dominate at T ~ 5.5 K. (d) R xx versus H plots in xy plane at different temperatures. 
The red lines are fits to equation (1). 
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FIG. 4: (Color online) MC of additional three Mn-Bi2Se3 samples in perpendicular field (H \\ c 
axis), (a) MC hysteresis is observed in the most lightly doped (Bi/Mn = 23.6) sample. Hysteretic 
MC jumps are observed in the heavily doped samples (Bi/Mn ^ 10.3), shown in panel (b) and (c). 
The arrows indicate the field sweep directions and the curves are shifted with respect to the curve 
measured at the highest temperature for clarity. 
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FIG. 5: (Color online) Hall conductivity of samples Bi/Mn = 23.6 and 12.5 at 0.5 K. The arrows 
indicate the field sweep directions. 
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FIG. 6: (Color online) MC of Mn-Bi2Se3 in field perpendicular to the sample plane, (a) MC of 
the most lightly doped sample (Bi/Mn = 23.6) shows a cross-over from weak anti-localization to 
weak localization at T ~ 3.0 K, consistent with the Tq obtained from SQUID measurements, (b) 
MC of the Bi/Mn = 12.5 sample shows larger positive MC, which survies up to ~ 5 K. (c) MC of 
the Bi/Mn = 10.3 sample, (d) MC of the most highly doped sample (Bi/Mn = 8.3). The insets 
in panel (a) - (d) plot temperature dependent Aa xx at fixed perpendicular magnetic field, (e) MC 
of four Mn-Bi2Se3 and one undoped Bi2Se3 sample, (f) Simulation using euqation (2) showing 
cross-over from weak anti- localization to weak localization as parameter A/2Ep changes. 
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FIG. 7: (Color online) (a) Core level spectroscopies on Mn-Bi2Se3 thin film before and after 
the decapping procedure, (b) High resolution ARPES dispersion mapping and corresponding 
energy dispersion curves (EDCs) of an undoped Bi2Se3 thin film along high symmetry M — V — M 
direction. The EDC that crosses the T-point is highlighted in green color, (c) High resolution 
ARPES dispersion mapping and corresponding EDCs of Bi/Mn = 15 Mn-Bi2Se3 thin film along 
high symmetry M — V — M direction. The EDC that crosses the T-point is highlighted in green. 
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